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00 I We investigate the dynamics of an initially disentangled Gaussian state on a general finite 

0^ ' symmetric graph. As concrete examples we obtain properties of this dynamics on mean field graphs 

^-P , (also called fully connected or complete graphs) of arbitrary sizes. In the same way that chains 

can be used for transmitting entanglement by their natural dynamics, these graphs can be used 

^5 ■ fo store entanglement. We also consider two kinds of regular polyhedron which show interesting 

f^ ' features of entanglement sharing. 
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1 Introduction 

One of the basic problems of quantum information processing is the problem of entanglement transfor- 
mation or more generally manipulation of entanglement. For long distances one usually uses photons 
through optical fibres or free air to transmit entanglement . However for short distances other methods 
are being explored which are based on using discrete sets of interacting quantum systems such as spins 
[D El O m [5] or harmonic oscillators [6j , which by their natural dynamics can generate and transmit 
entanglement. In particular in [6], one dimensional lattices of harmonic oscillators coupled by two 
different types of Hamiltonians, were studied and various phenomena were investigated with regard to 
entanglement generation and transmission. Among other things it was shown that the largest amount 
of entanglement between two oscillators is always obtained when one places them at the two ends of 
an open chain. This maximality was attributed to the fact that in this case the two oscillators have 
fewer neighbors to which they can become entangled. It was argued in [6j that besides linear arrays 
of oscillators, other geometries, in principle any arrangement corresponding to weighted graphs are 
worth of study, since they can act as building blocks of more complicated networks. In [6] itself two 
other geometries, namely a Y shape geometry which mimics a beam splitter and another geometry 
corresponding to an interferometer were studied. 

In this article we want to extend these considerations in one particular direction, namely we want 
to study compact and symmetric geometries, i.e. symmetric graphs of finite size. The basic motivation 
is that in contrast to the geometries considered in [71 E] which were suitable for transmission of en- 
tanglement, finite graphs are suitable for storing entanglement. As any other resource, entanglement 
needs to be stored for use in later suitable times and hence in any complicated network, building 
blocks which can store entanglement, should be implemented. In the simplest electrical analogy we 
may think of finite geometries as capacitors and linear arrays of the type considered in [6] as resistors 
or transmission lines. 

However in contrast to the static properties of entanglement, for which various symmetric graphs have 
been considered [8], for our purpose, only one type of symmetric graph seems to be useful, namely the 
mean field or a fully connected graph. The reason is the very simple temporal behavior of entangle- 
ment on these graphs, compared with the complicated behavior of arbitrary graphs. In fact a system 
of harmonic oscillators on a mean field graphs has only two natural frequencies which makes the result- 
ing time development of entanglement quite simple and easily controllable, while for other symmetric 
graphs, this is not the case. If as in [1] we are to extract entanglement at an optimal time, then it is of 
utmost importance that the dynamics of entanglement follows a simple and not a complicated pattern. 

For that reason we mostly consider mean field clusters of arbitrary sizes and determine how an 
originally disentangled set of harmonic oscillators positioned on the nodes of such a cluster, when 
coupled to each other, develop a pairwise entanglement between themselves. How this entanglement 
develops in time, what is its maximum value, and how it depends on the size of the cluster. We stress 
that our general setting is apt for analysis of any symmetric graph and we indeed include two other 
graphs for observing some other phenomena. 

The structure of this paper is as follows: In section [2] we briefly review the Gaussian states and 
their entanglement properties, especially we remind the closed formula for entanglement of Formation 
(EoF) [9] which in the context of symmetric graphs is more suitable than negativity as a measure 



of entanglement. In section [3] we study the dynamics of a Gaussian state on an arbitrary symmetric 
graph and obtain closed formulas for the EoF between any two sites as a function of time. This for- 
mula reduces the calculation of the EoF to the diagonalization of the adjacency matrix of the graph. 
In section m we study in detail the simplest graph consisting of two nodes. In section [5] we specialize 
to the mean field graphs where our concrete results are reported in figures ([51 [21 HD and table (0) . 

Finally in section ^ we consider two other symmetric graphs for comparison and draw some 
conclusions about the sharing of entanglement which challenge the arguments of [6] on this issue. 

2 Preliminaries on Gaussian States 

In this section we collect the rudimentary material on Gaussian states that we need in the sequel. 
References [lOl [11] can be consulted for rather detailed reviews on the subject of Gaussian states. 
Let R := (xi,a;2, • • • X]\f,Pi,P2, • • -Pn) be A^ conjugate operators characterizing N modes and subject 
to the canonical commutation relations 

[Rk,Rl] = i(^kh 

where o" = I j "is the 2n dimensional symplectic matrix and I„ denotes the n dimensional 

unit matrix. 

A quantum state p is called Gaussian if its characteristic function defined as C(^) := tr(e^*^*'^''' 'p), 

is a Gaussian function of the ^ variables, namely if 

where we have assumed that linear terms have been removed by suitable unitary transformations. The 
matrix F, called the covariance matrix of the state, encodes all the correlations in the form 

Tki := {RkRi + RiRk) - 2{Rk){Ri). 

For a two mode symmetric Gaussian state, one in which there is a symmetry with respect to the 
interchange of the two modes, the covariance matrix will be 

(1) 

where the modes have been arranged in the order xi, pi, X2, P2 and a and (3 are 2x2 symmetric 
matrices. By symplectic transformations the covariance matrix of a two mode symmetric Gaussian 
state can always be put into the standard form (in the order xi,X2,pi,P2) 




Ti rCp 

V kp n J 



(2) 



where k^ > > kp and kx > \kp\. The entries of the standard form of Tg can be determined from the 
fohowing symplectic invariants: 

n^ = det a (3) 

kxkp = det/3 (4) 

{n^ - kl){n^ - kl) = detT. (5) 

For a symmetric Gaussian state a closed formula for the entanglement of Formation has been derived 
in [9]. Note that there are other criteria for studying the entanglement or separability of Gaussian 
states [H], however we use only Entanglement of Formation here to take advantage of the inherent 
built-in symmetry of our graphs. EoF of a Gaussian state p, denoted simply by E{p) is expressed as 
follows: 

E{p) := C+ log2 C+ - C_ log2 C_ , (6) 

in which 

r (1^^)' (7^ 

and 

A := min{l, 6 := \/{n — kx){n + kp)). (8) 

Thus a state is entangled only \i 5 < 1. 

Note that 5 can be expressed in terms of the original covariance matrix. To express it we denote 



u 

V 

w 



det a (9) 

det/3 (10) 

detr (11) 



and 

Then a simple calculation gives 



^■-U^ +y^ -yj, (12) 



,-^ = „_„_^i-^lEi»_^i±^^liE15!. (13) 

In the following sections we use this equation for calculating the entanglement of a Gaussian state 
which is initially disentangled and evolves in time under a quadratic hamiltonian. 

3 Dynamics of entanglement of Gaussian states on symmetric graphs 

Consider a symmetric graph, having n-nodes, corresponding to an adjacency matrix A and a system 
of bosonic modes corresponding to the vertices of this graph interacting by a quadratic Hamiltonian. 
In this paper we consider a Hamiltonian of the form 

1 " 1 

'■=1 ihj) 



where the sum runs over adjacent nodes and c is a couphng constant. This Hamiltonian describes 
a simple mass-spring system of the form first studied by Plenio [6] in the context of entanglement 
dynamics. The above Hamiltonian can be written in the compact form 

H-ln^(y,)n (.5, 

where V and T (here equal to I) are the potential and the kinetic matrices. We include the case 
of arbitrary T (but commuting with V) for generality, since some other Hamiltonians like the one in 
[7] can be expressed in this way. 

The dynamics of R is easily determined by solving the equations of motion 

^ = -^[R,H] (16) 



or using [Ri, Rj] = ia^ 



«j 



dt \ -V 

The solution of the above equation is given by 



"" .MR J », Dr. (17) 



R(t) = expd "j^ ^ jt)fl(O). (18) 

The explicit form of the evolution matrix is found by writing it as 

U{t) := exp{ ( ^^ l)*) = e5-^®(^-^)+f-.«(^+^). (19) 

In order to find the explicit form of the evolution matrix we use the following 
Lemma: For any two commuting matrices A and B, the following identity holds: 

^Am+-B^a. = J ^ cos VAB + I a+ W 4 - fT_ W ^ I sin VAB. (20) 

This lemma is proved by a simple application of the identity e'^"'^ = cos 9 + i sin Ofi • a. 
Using the above lemma, we find the final form of the evolution matrix 



U{t) = I cos VVTt + ( ^/o-+ — - o-_ J— I sin VVT 




(21) 

For the case we consider in this article the kinetic matrix T is identity (T = /) and so with the 
definition W := VV, U{t) simplifies to 



TT^,\ T Txr, f Tir-1 TiA • Tir, / COSWt W ^ SmWt \ ,„„, 

U{t)=I(E)cosWt+[a+®W ^ -a-®W]smWt=\ ^_. . ^_ ^_ . 22 

^ ' \ J \ —WsmWt cosWt I ^ ' 

From the definition of the covariance matrix we find 

r(t) = t/(t)r(o)[/^(t). (23) 

Let us consider the case where the initial state is a completely uncorrelated state with r(0) = /. 
The covariance matrix as a function of time will then be given by 

vu\ - nmuTm _ ( c°«' ^* + ^"' ^i^' ^* (^'' " ^) ^in Wt cos Wt \ 

n^)-^ {^) ^ \.^)-y (^-1 _ ^) g- j^ ^^ ^Qg ^^ ^Qg2 ^^ ^ ^2 g- j^ ^^ j • {^V 

The covariance matrix between any two modes (sites of the graph) is determined by extracting only 
the sub-matrix pertaining to those two sites. For this we need the matrix which diagonalizes V. Let 
Q,~^WU = Wd where Wd is a diagonal matrices with diagonal elements cjj. The matrix Q. is the 
matrix which diagonalizes the adjacency matrix of the graph. 



r(t) = {^®n){ % J!| ) (^ e ^f (25) 



Yu pi 

xp pp 



where 



^xx = diag (cos u>it + uj^ sin ujit) 
Tpp = diag (cos Wjt + uj^ sin uJit) 

F = diag {{co^^ — uJi) sin ujit cos uJit). (26) 

Then the covariance matrix between any two modes say modes 1 and 2 will be the form 



ra;i,a;i 


— ra;2,X2 




ra;i,a;2 


= ra;2,a;i 


n 
= ^^li^2j(ra;x)i 


r 


— r 


n 

= ^J7ii^(F^)i 


''Pl,P2 


— J- P2,Pl 


n 



1=1 



J-aiijpi — ^ X2,P2 ~ 2^ ^i^ ■ 



xp) 
i=l 



-*- a;i,p2 -*- a;2,Pi / , "lt"2tl.-*- xv)i- l^'J 

i=l 

Therefore in each case we should only determine the matrix which diagonalizes the potential 
matrix V and from (j26l and [27| ) determine the eigenvalues Wj . In the forthcoming sections we use 
this formalism to determine the dynamics of entanglement between any two modes on a wide variety 
of symmetric graphs. Note that the entanglement of formation is defined only for symmetric Gaussian 
states, and in this paper we are considering only symmetric graphs. So in all of the graphs that we 
consider, this entanglement is invariant under isomorphism of graphs. 



4 The simplest example, A Two-Mode System 

As the simplest example we consider a two mode system represented by a simple graph consisting of 
two nodes and a link connecting them. 
The Hamiltonian is 

^ = \{v\+vi) + \iA + ^i + c(xi - X2)2), (28) 



which corresponds to the matrices 



aniV=[ 1 _+; J-;; . (29) 



The eigenvalues of the matrix V are readily obtained to be 1 and uj := \/l + 2c. 
The covariance matrix is given by 



" = 7l(5 -i)- P») 



and 



^xx — diag{l, cos'^ ujt + uj^^ sin^ ujt) 

r = diag{l, cos'^ ujt + u)'^ siu^ cot) 

^x,p — diag{0,{uj^ — uj) sin oot cos u!t). (31) 

Using (J13p we find the parameter 5 which is essential for calculating the entanglement of the state. 
The result is 

1 + jiuJ-'^ - uj)"^ sin^ cot I (a;"^ - lo) sinujt \ . (32) 

The entanglement of the state is shown in figure ([T]) for two difference coupling constants. 
The entanglement oscillates at the natural frequency to. Its maximum value is achieved for the 
minimum value of 6, or 



5{ut = ^) = Vl + \(^'' - ^)' - ^ I (^"' - ^) 




t 

Figure 1: (color online) The entanglement between the two modes in a simple two- vertex graph, as a 
function of time, for two different couplings, c = 2 (dashed line) and c = 8.3 (solid line). 

As we increase the frequency or the coupling constant, the entanglement becomes flat in most of the 
period and develops cusp singularities in half periods. Note that 6 ranges between (for lo — > cxd) and 
1 (for LO = 1). Thus the maximum entanglement increases unboundedly by increasing the strength of 
the interaction c. In fact the maximum entanglement increases as log cu for large coupling constants 
c. 

Thus a two- vertex graph can be used as a storage device for entanglement the "capacity" of which 
increases logarithmically with the coupling constant c. Moreover as the flatness of the curve in figure 
([T]) shows, for very large coupling constants we can extract this maximum entanglement at any time 
we wish except for a discrete set of points. 



5 Mean Field Clusters 

We now consider a mean field cluster of N vertices in which every vertex is connected to A^ — 1 other 
vertices. The adjacency matrix for a mean field graph is given by 



A = E-L 



(33) 



where E is the matrix all of whose entries are equal to 1, Eij = 1 V i, and j. The potential 
matrix of this graph is given by 

V = {1 + Nc)I - cE. (34) 



The matrix E and hence V can easily be diagonalized. We have 

Eeo = Neo, Eek = 0, k = 1,2,- ■ ■ N - I, 



(35) 




t 

Figure 2: The re-scaled entanglement as a function of time for a mean field cluster of size A^ = 20 at 
c = 1, the time spans one full period. 



where 



eo 
efc 



(l,l,---lf, 



1 



= (1,1,1, 



VHk + T) 

Thus the eigenvalues of V will be given by 

Wo = 1 



k,---o) 



T 



k = l,2,---N -1. 



UJl 



■LON-l 



ViT1^=:lo. 



(36) 



(37) 



The eigenvectors eo to ctv-i derived above easily yield the diagonlizing matrix 0, {0,ij = {ej)i) from 
which we can obtain after straightforward calculations from (j26p and ()27p the following parameters of 
the covariance matrix between any two sites say sites 1 and 2: 



r 

xi,a;i 

r 

^ P1,P2 
■>- xi,pi 
■■■ Xl,p2 



1 



1 



^ X2, X2 = ir + i'^ - jr) (cos^ ^t + Lo ^sin^wt) 

r^a.xi = — (1 - w"^) sin^ wt 

1 1 

rp2,P2 = — + (1-— )(cos2u;t + a;2sin2wt) 

^P2,Pi = ^(l-w^)sin2a;t 

rj^2,P2 = (1 — at'^^^' — u>) sin cut cos ujt 

Tx2,pi = ~ir^{^~ — uj) sin ut cos ujt 



(38) 




N 

Figure 3: The maximum re-scaled EoF for mean field clusters for different coupling constants. 



One can obtain the standard form of this matrix by using the symplectic invariants. They read in 
the present case 



u 



w 



i^c '"• •"'"-^"" 



A^ 



)(u; — u; ) sin oot, 



-1 



{lo — LO ) sin ujt, 



1 + 2{lo -u)' 



^l\2 



.N 



Ar2 



I sin ujt. 



(39) 
(40) 
(41) 



Inserting these values in 1131 gives the entanglement for these graphs. Following [ISJ, we define 
re-scaled entanglement {E^-) which is A^ — 1 times the entanglement between any two nodes. This 
definition stems from the fact that a node shares its entanglement with its neighbors which are A^ — 1 
in number. 

Figure ^ shows the re-scaled entanglement for a mean field cluster of size 20 as a function of time. 
Figure ([3]) show the maximum re-scaled entanglement for mean field clusters as a function of their 
size for different coupling constants. 

For any cluster of size N and coupling constant c, the entanglement oscillates at a frequency 
oj := 2\/l -|- Nc and most of the time the two modes have an appreciable amount of entanglement. 
With increasing the coupling constant c, the amplitude of oscillation increases and saturates to a 
finite value for very large c, as long as A^ > 2 (figure H]) . Table ([5]) shows this saturated amplitude for 
clusters of different sizes. 



Maximum E 




Figure 4: The amplitude of oscillation of entanglement (E) in units of 0.01 ebits, for mean field clusters 
of different sizes as a function of coupling constant. The saturation values are always less than the 
bounds found in \7]. Note that for all the sizes A^ > 2, as c increases, E saturates to finite values. 

6 Sharing of entanglement 

As mentioned in the introduction, in [6] it was shown that in a linear lattice, the largest amount of 
entanglement between pairs of sites with the same distance, occurs when these two sites are at the 
end points of the lattice. This was attributed to the fact that the endpoints of the lattice have fewer 
neighbors to which they share their entanglement. In this regard it is instructive to consider two 
symmetric graphs which are not fully connected. The two graphs which we study are a six- vertex 
graph in the shape of octahedron and an eight-vertex graph in the shape of a cube. They are shown 
in figure ([5]) with numbered vertices. 

In the octahedron there are essentially two types of pairs, represented by the pair (1,2) and the pair 
(1,6). In each pair the number of neighbors of each node is the same. However the pair (1,6) although 
more apart than the pair (1,2) develops a much higher entanglement, figure[6l In terms of the number 
of edges, the distance between the nodes 1 and 2 is one, and there is only one shortest path which 
connects these two nodes, while the distance between nodes 1 and 6 is two, however there are four 
such shortest paths which connect these two nodes. Therefore it seems that entanglement between 
two site is not only affected by the number of their neighbors, but also by the number of shortest 
paths which connects these two sites to each other. To test this idea, we study the cube, which has 
three types of pairs represented by the (1,2), (1,6) and (1,7), with distances respectively given by 1, 2 
and 3 and the number of shortest paths respectively given by 1, 2, and 6. The entanglement is shown 
in figure ([7]) which confirms our assertion. Here we see a competition between the two factors. If we 
interpret the entanglement as a direct measure of quantum correlations then these two figures show 
a very intriguing property of entanglement: there are times where remote sites are strongly quantum 
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Maximum Er 


2 


CO 


3 


0.803 


4 


0.592 


5 


0.484 


6 


0.415 


7 


0.365 


8 


0.328 


9 


0.298 


10 


0.274 


15 


0.196 


20 


0.156 


30 


0.113 



Table 1: The saturated amplitude of re-scaled entanglement for mean field graphs of different sizes. 

1 





Figure 5: (color online) Two symmetric graphs, Cube and the Octahedron, the entanglement between 
different pairs of vertices are plotted in figure figure [6] and [71 

correlated while the nearest sites have a small quantum correlation. 
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Figure 7: (color online) The pairwise entanglement between the sites (1,2), (1,3) and (1,6) in the 
Cube as a function of time. 
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